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CMB Maps at 0?5 Resolution 
I: Full-Sky Simulations and Basic Results 

G. Hinshaw 1 ' 2 , C.L. Bennett 3 , and A. Kogut 1 
ABSTRACT 

We have simulated full-sky maps of the Cosmic Microwave Background (CMB) 
anisotropy expected from Cold Dark Matter (CDM) models at 0?5 and 1?0 angular 
resolution. Statistical properties of the maps are presented as a function of sky 
coverage, angular resolution, and instrument noise, and the implications of these results 
for observability of the Doppler peak are discussed. The rms fluctuations in a map 
are not a particularly robust probe of the existence of a Doppler peak, however, a full 
correlation analysis can provide reasonable sensitivity. We find that sensitivity to the 
Doppler peak depends primarily on the fraction of sky covered, and only secondarily 
on the angular resolution and noise level. Color plates and one- dimensional scans of 
the maps are presented to visually illustrate the anisotropies. 

Subject headings: cosmic microwave background — cosmology: observations 

1. Introduction 

Observations of the cosmic microwave background anisotropy are fundamental to 
understanding the formation and evolution of structure in the universe. On angular scales 
larger than a few degrees, the CMB anisotropy traces the primordial density distribution, while 
measurements at higher resolution (9 <J 1°) probe the physical scales and causal mechanisms 
responsible for currently observed structure, offering greater power to distinguish between 
competing models of structure formation (Efstathiou, Bond, & White 1992; Coulson et al. 1994) 
Several groups have reported detections of anisotropy on degree angular scales (Schuster et al. 
1993; Wollack et al. 1994; Cheng et al. 1994; De Bernardis et al. 1994; Dragovan et al. 1994; 
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Devlin et al. 1994; Clapp et al. 1994) but the implications of these experiments for models of 
structure formation are unclear. 

Interpretation of observational results in a specific patch of sky are complicated by the fact 
that cosmological models predict the CMB anisotropy to be a single realization of a stochastic 
process whose properties (power spectrum, phase correlations, etc.) are predicted only for an 
ensemble average of equivalent realizations. Any single realization may be expected to vary 
about the ensemble average ("cosmic variance"). The generally small sky coverage of most 
existing observations (typically 10~ 3 to 10~ 4 of the full sky) provides another complicating 
factor: observations of small patches of the sky may not be representative of the sky as a whole, 
particularly in non- Gaussian models with distinctive rare features ("sample variance"). 

Although a large body of work exists providing analytic or Monte Carlo derivation of ensemble 
average properties of various models, relatively little has been published on the phenomenological 
properties expected for single medium-scale experiments. Questions of importance both to 
observers planning future instruments and to theorists analyzing existing data include the 
following: what is the role of sample variance in statistical analyses of small regions of the sky? 
What fraction of the sky must be observed to constrain models at specified confidence, particularly 
against Type II errors, accepting a hypothesis as true when it is false? What is the relation 
between limiting sources of uncertainty: the signal-to-noise ratio per pixel, cosmic variance, sample 
variance, and sky coverage? Are large-area scans with relatively greater instrument noise per pixel 
better than small-area deep observations? To what extent do features of the standard model such 
as the Doppler peak appear to observers as anomalies like isolated point sources? How important 
is it to produce a map of the microwave sky, as opposed to observing chopped scans? 

Some of these questions are beginning to be addressed (Scott, Srednicki, & White 1994; Luo 
1994; White 1994; Bond 1994) but many remain unanswered. In this Letter, we describe Monte 
Carlo simulations of CMB anisotropy in CDM models at 0?5 and 1?0 angular resolution, which 
we use to quantify the power of simple statistical tests that can be applied to the maps. As an 
example, we construct and evaluate a likelihood function for the universal baryon density, and 
show that the sensitivity to is determined primarily by the fraction of sky covered, and, to a 
lesser extent, by the angular resolution, and the instrument sensitivity. Further results on the 
analysis of unresolved features are reported in a companion paper (Kogut et al. 1995). 

2. Map Simulations 

We have simulated full-sky CMB maps at 0?5 and 1?0 angular resolution, Gaussian full width 
at half-maximum (FWHM), for a range of cosmologically interesting CDM parameters. The 
anisotropies expected in such models are expressed in terms of the mean angular power spectrum, 
specified by multipole amplitudes C'i, which, for Gaussian fluctuations, completely determine the 
enesemble- average properties of the sky. Figure 1, after Stompor (1994), shows the power spectra 
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of the CDM models we consider; in all cases we adopt a scale-invariant primordial spectrum of 
fluctuations (n = 1) with a quadrupole normalization of Q rms -PS = 20 /uK based on the two-year 
COBE-DMR data (Bennett et al. 1994; Gorski et al. 1994; Wright et al. 1994; Banday et al. 
1994). The dashed lines in the Figure represent the spectra convolved with 0?5 and 1?0 beams. A 
familiar feature in these spectra is the existence of the so-called Doppler peak at I ~ 220. For a 
given primordial Sachs- Wolfe spectrum and a given Hubble constant, the amplitude of the first 
Doppler maximum increases with the baryon density In this Letter, we restrict attention to 
flat models with S7o = 1, but it should be noted that the location of the first Doppler maximum, 
which probes the angular size of the horizon at the epoch of last scattering, depends primarily 

1/2 

on Oo : t max ~ 200 / O (Kamionkowski, Spergal & Sugiyama 1994). Thus, while a 1° beam is 
adequate for probing the Doppler peak region in a flat universe (see §4.), it is marginal in a open 
universe. 

To generate a particular realization of the CMB temperature observed at an 
instrumental dispersion at, we evaluate the sum over spherical harmonics, AT(0,(f>) = 
J2im a lmYi m (0,(f))e~2 i ( i+1 } a b , where the amplitudes ai m (2 < £ < 500) are Gaussian random 
variables of zero mean and variance C'i, ({(i£m a £i m i) = CiSfjj S mm i) appropriate to the model under 
consideration (Stompor 1994; Bond & Efstathiou 1987). Each map is pixelized according to the 
quadrilateralized spherical cube projection (White & Stemwedel 1992), which projects the entire 
sky onto six cube faces. Each face is divided into 2 2 ( w_1 ) approximately equal area, square pixels 
where (N — 1) is the index level of the pixelization. This construction allows for very efficient 
analysis of, for example, subsets of high resolution pixels contained within a given low resolution 
pixel. The current maps are pixelized with N = 9, which gives 393216 pixels of size ~ 0?32 X 0?32. 

Figure 2 shows color images of selected maps at 0?5 angular resolution to illustrate visual 
features of the anisotropies expected in CDM models. Each map was generated with an 
independent random number seed, so the large scale features vary from map to map, but in 
all cases the structure due to the lowest order multipoles (I ^ 10) is readily apparent, and 
statistically equivalent. More striking differences are apparent on smaller angular scales where the 
power spectrum exhibits significant structure. The upper map was generated with a flat spectrum, 
defined as C'l/iir = 1.2 Q 2 rm s-Ps/W^ + !))> while the middle and lower maps were generated with 
fi& = 0.05,0.20 CDM spectra respectively. The increase in high frequency power evident from top 
to bottom is a plain manifestation of the increase in Doppler peak amplitude. In fact, the high 
frequency power in these maps is, in some cases, visually remniscent of a population of unresolved 
sources in the microwave sky. We explore this subject in a companion paper (Kogut et al. 1995). 
The peak amplitudes in these 0?5 resolution maps typically range from ~ 200 /uK in the flat 
spectrum model to ~ 400 in the fi 6 = 0.20 CDM model. 

Most present medium-scale experiments present their data in the form of one-dimensional, 
differential scans of a strip of sky. We crudely mimic such scans in Figure 3 where the left-hand 
panels show selected map profiles along a 10° strip of sky near the north pole, and the right-hand 
panels show the corresponding single difference observations using a beam separation of 1?28. (We 
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have not subtracted a best fit gradient from the difference scans, as many observers do.) Note 
that relatively large and seemingly isolated features are not uncommon in these profiles. 

3. rms Analysis 

A basic measure of anisotropy in the maps is the pixel-pixel rms temperature fluctuation. 
Table f gives the rms fluctuations observed in the simulated maps as a function of the baryon 
content of the universe, the solid angle of the region selected for observation, and the angular 
resolution of the observations. The first row in each table gives the expected rms, and its standard 
deviation, for a full-sky map generated from the corresponding model parameters. The second 
row lists the rms observed in each individual full-sky map. In no case does the deviation from the 
ensemble average exceed 2a. The third row lists the mean rms observed in six patches of the sky, 
each subtending |th of the sky, as defined by the six pixelization cube faces. The uncertainty in 
this row is the standard deviation of the six observations and is given to illustrate the typical level 
of sample variance that can be expected in such observations. Successive rows in the table list the 
mean rms fluctuations observed in successively smaller (and more numerous) portions of the sky. 
In particular, the final row presents results for the 1536 patches of size ~ 5° X 5° (8.2 X 10~ 3 sr), 
which is comparable to, or larger than, the coverage in most current medium-scale observations. 

The trend towards decreasing rms values with decreasing solid angle merely reflects the fact 
that small patches are insensitive to features whose angular scale is larger than the patch size. 
Such structure is subtracted off as part of the mean temperature in the patch. The increasing 
scatter with decreasing solid angle reflects the corresponding growth of sample variance. Note 
that the smallest patches considered here have a sample variance on the order of 15 — 20%, which 
is comparable to the cosmic variance in the full-sky COBE maps at low I. The trend towards 
increasing rms with increasing is a result of the growth of the Doppler peak with The 
scatter in the observed rms provides a crude means for estimating how well the height of the 
Doppler peak can be resolved for a given angular resolution and sky coverage, in the absence of 
experimental uncertainty. For example, a map covering |th of the sky with 1° resolution (see 
Table 1) can resolve a Doppler peak characterized by = 0.05 from a model with no Doppler 
peak, but a map covering 20° X 20° cannot, based on rms values alone. The situation is much 
better at 0?5 resolution where, even in the smallest patches, the difference between an = 0.05 
CDM sky and a flat spectrum sky is roughly 3(7. 



4. Likelihood Analysis 
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By definition, the rms fluctuations in a map are only sensitive to the pixel variance which, in 
turn, depends only on the total power within the window function of the experiment 



where t{ is the temperature in pixel i of a map, the angled brackets denote a universal ensemble 
average, Bj = e~^ +1 ) <7 f> is the beam factor and Ci is the power spectrum. However, the angular 
power spectrum also predicts the correlations expected between pairs of pixels and, for Gaussian 
fluctuations, the covariance matrix 
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fully characterizes the map. Here P;(cosa 8 'j) is the Legendre polynomial of order £, and aij is the 
angle between pixels i and j. The likelihood of observing a map with pixel temperatures /, given 
a power spectrum C;, is 

where N is the number of pixels in the map. Assuming a uniform prior distribution of model 
parameters, the likelihood of a model, given a map /, is then 

L?r.M-HC e yt 

L(Ce) oc — . 
V Vdet M(C t ) 

In principal, the additional information contained in the likelihood function should provide 
superior sensitivity to the such features as the height of the Doppler peak, as parameterized by fife. 
To test this, we evaluate the above likelihood as a function of for the two smallest patch sizes 
considered above at 0?5 and 1?0 resolution. More precisely, we simulate 100 independent patches 
of size 10?4 X 10?4 and 5?2 X 5?2 at each resolution using the parameters Q rms -ps = 20 /iK, 
n = 1, H = 50 km s" 1 Mpc" 1 , and tt b = 0.05. Each 10?4 (5?2) patch contains 1024 (256) square 
pixels. We simulate uncorrelated receiver noise in the maps by adding a random Gaussian number 
to each pixel with a fixed standard deviation (this implies uniform sky coverage in the patch, we 
then modify the above covariance matrix to account for this contribution: Mij — ► Mij + oAj, 
where a is the noise standard deviation per pixel). According to Table 1 the rms fluctuations due 
to an fife = 0.05 model are typically 50 — 70 [iK for these patches, so we add receiver noise ranging 
from 5 — 40 [iK per pixel to probe a noise-to-signal regime from ~ to ~ 1. 

Figure 4 shows the likelihood functions obtained from single, randomly selected patches of 
the sky for each resolution and patch size. The width of the likelihood function, which measures 
the sensitivity to is nearly independent of the angular resolution of the map (in the small 
range considered), but strongly sensitive to the fraction of sky observed. This point is reiterated 
in Table 2 which gives the statistics of the recovered maximum likelihood values for the ensemble 
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of patches. In each case, the median of the 100 determinations is always within 2a of the input, 
fi& = 0.05, while the standard deviation ranges from 0.12 in the large, high-resolution, low-noise 
patches, to 0.41 in the small, high-noise, low-resolution patches. This table also demonstrates 
that in the regime of moderately high signal-to-noise, a <J 40 /iK, the sensitivity to is nearly 
independent of noise level. We interpret these results as follows: a precise determination of the 
height of the Doppler peak requires a large number of "typical" high spatial frequency (I ~ 200) 
peaks and troughs to determine a "typical" amplitude. Moreover, given sufficient instrument 
sensitivity, a 1° beam will not dilute these features to the point of non-detection. We have also 
evaluated the above likelihood function, in each of our simulated = 0.05 patches, for the flat 
spectrum model. In every patch the likelihood of this model is exponentially suppressed relative 
to the standard CDM models. 

Of course the above analysis makes the optimistic assumption that all other cosmological 
parameters are known exactly and that there are no significant systematic effects (see, e.g., 
Wilkinson 1994) or foreground contaminants (see, e.g. Bennett et al. 1992). There is also 
considerable degeneracy amongst cosmological parameters with regard to the interpretation of the 
CMB anisotropy power spectrum in the sense that different parameter combinations can yield 
similar angular power spectra (Bond et al. 1994). The conclusions we draw about experimental 
sensitivities to are subject to such considerations, however our conclusions regarding the 
experimental sensitivities to the CMB power spectrum itself are valid independent of theoretical 
models. Our general conclusion is that, given limited observing time, optimal sensitivity to the 
CMB power spectrum is obtained by covering as large a portion of the sky as possible consistent 
with obtaining a moderately high signal to noise ratio. 

We are very grateful to Radek Stompor for providing us with CDM power spectra in advance 
of publication. Also, it is a pleasure to acknowledge useful conversations with Tony Banday and 
Krys Gorski. 
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Table 1. pixel-pixel rms fluctuations (/iK) 



Patch size a 


flat c 


n b = 0.01 


Model 
tt b = 0.05 


fi 6 = 0.20 






0?5 resolution 






ensemble d 


66.7 ± 2.5 


80.4 ± 2.1 


85.7 ± 2.0 


105.9 ± 1.6 


full sky 


62.3 


82.5 


88.1 


105.9 


~ 90° X 90° 


62.0 ±2.7 


81.4 ± 5.7 


86.0 ±4.3 


105.7 ± 2.3 


~ 45° X 45° 


60.2 ±5.7 


77.2 ±5.1 


82.2 ±4.1 


104.8 ±3.6 


~ 20° X 20° 


55.7 ± 6.4 


72.8 ±5.6 


78.9 ± 5.4 


100.7 ±4.8 


~ 10° X 10° 


49.9 ± 7.3 


67.6 ±6.9 


73.9 ± 6.5 


96.5 ± 6.4 


~ 5° x 5° 


42.8 ±8.0 


61.4 ± 8.2 


68.5 ± 9.1 


92.1 ± 10.1 






1?0 resolution 






ensemble 


61.5 ±2.7 


68.2 ±2.5 


71.2 ± 2.3 


81.0 ± 2.1 


full sky 


58.8 


70.3 


72.5 


80.7 


~ 90° X 90° 


58.3 ±6.7 


68.5 ±5.0 


72.1 ± 3.3 


78.7 ± 2.3 


~ 45° X 45° 


55.8 ± 7.0 


63.3 ±4.2 


66.8 ± 6.0 


76.9 ± 5.0 


~ 20° X 20° 


50.0 ± 7.3 


58.3 ±6.9 


61.2 ± 5.6 


73.5 ± 6.3 


~ 10° X 10° 


43.1 ± 7.9 


52.0 ± 7.1 


55.4 ± 7.3 


67.7 ± 7.5 


~ 5° x 5° 


35.1 ± 8.7 


44.6 ±8.5 


48.7 ± 9.0 


61.3 ± 10.7 


a The patch 


sizes listed above are 


approximate, for simplicity. 


The exact patch 


sizes, based on 



solid angles subtended by the low resolution pixels, are: 47r, 47r/6, 47r/24, 47r/96, 47r/384, 47r/1536. 

b All other model parameters for this Table were fixed: Q rm s-PS = 20 /iK, n = 1, and Ho = 50 
km s _1 Mpc -1 . 

c The flat model is defined as CejAir = 1.2 Q 2 rms _ PS /(£(£ + 1)). 

d The ensemble average rms and its standard deviation expected from a given model. 
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Table 2. Derived 


likelihood maxima 


from = 0.05 simulations 


a 


Patch size 


5 


10 


Noise per pixel (/«K) b 
20 


40 


5?2 x 5?2 
10?4 x 10?4 


0.051 ± 0.025 
0.048 ± 0.012 


0?5 resolution 
0.051 ± 0.027 
0.051 ± 0.014 


0.053 ±0.030 
0.051 ± 0.014 


0.055 ± 0.032 
0.048 ± 0.017 


5?2 x 5?2 
10?4 x 10?4 


0.046 ± 0.027 
0.049 ± 0.015 


1?0 resolution 
0.047 ± 0.030 
0.048 ±0.017 


0.049 ±0.037 
0.049 ±0.018 


0.042 ± 0.041 
0.046 ± 0.021 



a Table entries are the median and rms of 100 maximum likelihood determinations for nominal 
values of the other parameters (see §4.). 

b The modeled receiver noise is uncorrelated from pixel to pixel. Note that for a fixed observing 
duration, the noise per pixel in a 10?4 patch would be twice that in a 5?2 patch. 
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Fig. 1. — Spectra of the simulated CDM models, after Stompor (1994). All spectra have a scale- 
invariant primordial spectrum (n = 1), and have been normalized to a quadrupole anisotropy of 
Qrms-PS = 20 /uK. The flat model corresponds to a pure Harrison- Zel'dovich spectrum, and is 
included as representative of models without a Doppler peak. The CDM models plotted have 
a Hubble constant H = 50 km s" 1 Mpc" 1 , and baryon density fi 6 = 0.01,0.03,0.05,0.10,0.20 
increasing from bottom to top. The dashed lines in the top (bottom) panel correspond to spectra 
convolved with a 0?5 (1?0) FWHM Gaussian beam. 
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Fig. 2. — Simulated full-sky CMB maps at 0?5 resolution. The temperature range, as indicated by 
the color bar, is ±300 [iK. top) Map generated with a flat spectrum model (see text), middle) Map 
generated using a CDM spectrum with baryon density = 0.05 and nominal values for the other 
parameters (see text), bottom) Same as middle, except = 0.20. 
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Fig. 3. — Selected map profiles (left) and corresponding single difference scans, with a 1?28 beam 
throw (right), a) Profile from the flat spectrum model, and its corresponding single difference scan. 
b) Same as a) with an = 0.05 CDM model, c) Same as a) with an = 0.20 CDM model. 
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0.00 0.05 0.10 0.15 0.20 

Fig. 4. — Likelihood functions for the baryon density for selected sky patches. In each case the 
noise per pixel was 10 /iK, and the input value of was 0.05. Note that the width of the likelihood 
depends relatively strongly on the patch size, but not on the resolution. 
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